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Abstract 

Inflationary solutions are constructed in a specific five-dimensional model with boundaries mo- 
tivated by heterotic M-theory. We concentrate on the case where the vacuum energy is provided 
by potentials on those boundaries. It is pointed out that the presence of such potentials necessarily 
excites bulk Kaluza-Klein modes. We distinguish a linear and a non-linear regime for those modes. 
In the linear regime, inflation can be discussed in an effective four-dimensional theory in the con- 
ventional way. We lift a four— dimensional inflating solution up to five dimensions where it represents 
an inflating domain wall pair. This shows explicitly the inhomogeneity in the fifth dimension. Wc 
also demonstrate the existence of inflating solutions with unconventional properties in the non-linear 
regime. Specifically, we find solutions with and without an horizon between the two boundaries. 
These solutions have certain problems associated with the stability of the additional dimension and 
the persistence of initial excitations of the Kaluza-Klein modes. 



1 Introduction 



Two important theoretical developments, the advent of M-theory and the discovery of branes have 
recently stimulated new ideas in early universe cosmology. There has been considerable activity on 
various cosmological aspects of M-theory over the past two years ][l. | |24|] . The cosmology of Hof ava- 
Witten theory p5], [ 



27, pq|, which describes M-theory on the orbifold S /Z2, however, is much 



less studied so far 



30, 31, 32, 33], This theory describes the strong coupling dual of the Eg x E% 
heterotic string and is, therefore, of great importance for M-theory particle phenomenology. Clearly, 
this property makes it a very interesting starting point for cosmology as well. 

Various aspects of branes might be important in early universe cosmology, such as their ability to 
smooth out singularities ||, |6|, O, 21, 22] and their thermodynamical properties [1C, 12, [34], J35| . Most 
obviously, however, they play a role in the cosmology of particle physics models that have branes in their 
vacuum structure and, more specifically, that lead to low-energy theories arising from the worldvolumes 



of branes. Such models appear in the context of brane boxes [|36|, |37j, type I string theory |38|, |39|, |40j, |41[| 
and M-theory on S 1 /Z2 fl42"| , 43, [i4|, 45 1. A characteristic feature of many of those models is the 
possibility of one or more compact dimensions being large compared to the fundamental length scale 
of the theory. Such a situation can be described by a Kaluza-Klein theory with gravity and possibly 
other fields in the bulk coupled to a four-dimensional "brane-like" object with the standard model 
fields on its worldvolume. A wide spectrum of scales has been proposed for such models. These range 
from fundamental scales around the GUT scale with the energy scale associated with the additional 
dimensions being an order of magnitude or so smaller, to models with a fundamental scale of order a 
TeV with macroscopic additional dimensions pq| . It is clearly interesting to explore the cosmology of 
these models and, recently, some work in this direction |47], 21, [H], 33, Q has been done. 

In this paper, we would like to study the important issue of how inflation relates to these new 
theoretical ideas. For recent related work on inflation see [^, 4£, 5C, 2^, 24, mj. Rather than presenting 
a general scenario, we will concentrate on a specific model which incorporates the M-theory as well as 
the brane aspects. This model can be interpreted as a part of the five-dimensional effective action of 
M-theory on S /Z% ^, 52, 43, 53] obtained by reducing the 11-dimensional theory on a Calabi-Yau 
three-fold. The five-dimensional space of this theory has the structure M5 = S jZ-i X M.^ and contains 



(i) 



I (2) 



two four-dimensional orbifold fixed planes (or boundaries) M\ and M\ . It consists of gauged J\f = 1 
supergravity plus vector- and hypermultiplets in the bulk coupled to N = 1 theories with gauge and 
chiral multiplets on the orbifold fixed planes. The vacuum solution of this theory |^2| is a BPS double 
three-brane (domain wall) with the three-brane worldvolumes identified with the orbifold planes. Upon 
reduction to four dimensions on this vacuum solution, one arrives at an N = 1 supergravity theory which 
is a candidate for a realistic particle physics model from M-theory. The hidden and observable fields 
in this model arise from the "three-brane orbifold planes". The theory, therefore, allows us to study 
cosmology in a potentially realistic particle physics environment and provides a concrete realization of 



the general idea of "getting the standard model from a brane". Simple cosmological solutions of this 



theory have been found in ref. [31, 33]. In ref. [44, 45] non-perturbative vacua of heterotic M-theory 



containing five-branes have been constructed. In the five-dimensional effective theory, these five-branes 
appear as three-branes which, in addition to the two orbifold planes, are coupled to the bulk. It would 
clearly be interesting to study cosmological solutions of these more general theories with additional 
three-branes. In this paper, however, we restrict ourselves to the simple setting of two orbifold planes. 

For the application to cosmology, we have (consistently) truncated this theory to a minimal field 
content suitable for inflationary models. Specifically, in the bulk we have kept gravity and a scalar field 
4> (the volume modulus of the internal Calabi-Yau space) with a potential V of non-perturbative origin. 
In addition, on each orbifold plane we have also kept a scalar field fa with a potential V&. The theory 
is characterized by three scales, namely the fundamental scale M5 of the five-dimensional theory, the 
separation R of the orbifold planes and a scale v that sets the height of certain explicit potentials for 
the bulk field <f>. These potentials are responsible for the existence of the domain wall solution. 

Hence, we have a very simple setting with one additional dimension and one "candidate inflaton" 
with potential in each part of the space. In addition to the M-theory relation, we, clearly, also have a 
simple starting point to study inflation in the general context of models with large additional dimensions. 

The goal of this paper is not to construct explicit inflationary models by choosing specific potentials 
for the scalar fields in the theory. Rather, we are interested in how the specific structure of the theory, 
that is, the coupling of the five-dimensional bulk to four-dimensional boundary theories, effects inflation. 
We distinguish two different types of inflation, which we call bulk (or moduli) and boundary (or matter 
field) inflation. For bulk inflation the vacuum energy is predominantly provided by the bulk potential 
V, whereas for boundary inflation the boundary potentials V& dominate. In this paper, we concentrate 
on the second case of boundary inflation. This option is particularly interesting in that it directly relates 
to the presence of the characteristic boundary theories. Moreover, inflation from matter fields seems 
to be in better accord with current directions in four-dimensional inflationary model building pll] than 
modular inflation. 

Let us summarize our main results. One of the main themes of this paper is that energy density 
on the orbifold planes provides source terms localized on the fixed points in the additional dimension 
and, hence, excites bulk fields. In particular, this applies to vacuum energy on these planes as needed 
for boundary inflation. Our first conclusion is that boundary inflation is necessarily inhomogeneous 
in the additional dimension, or, in other words, excites Kaluza-Klein modes. The magnitude of those 
excitations is controlled by the dimensionless parameter £» = V^R/M^. For |ej| <C 1 the excitations 
can be described by linearized gravity. This approximation breaks down if |ej| ^> 1. One then has to 
use the full non-linear theory. Using the COBE normalization and the typical magnitude of heterotic 
M-theory scales, we argue that inflation in this theory may take place in both regimes. Interestingly, if 



the orbifold and Calabi-Yau scales during inflation are at their physical values [28], the theory becomes 
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linear for V±/ 4 <C 6 • 10 16 GeV, right below the COBE scale. 

The proportionality of £j to the size R of the additional dimension can be understood from the 
linear behaviour of the one-dimensional Green's function. For more than one additional dimension, 
the Green's function is logarithmic or follows an inverse power law. As a consequence, in those cases, 
the theory is always in the linear regime. The case of heterotic M-theory with one large dimension is, 
therefore, the only one where inflation may take place in the non-linear regime. 

In the linear regime, we can compute a sensible four-dimensional effective theory by integrating 
out the Kaluza-Klein modes induced by the boundary sources. We show explicitly how this leads to 
corrections in the four-dimensional theory. Our basic statement is that, in the linear regime, inflation 
can essentially be treated in the effective four-dimensional jV = 1 supergravity theory. Nevertheless, 
to get a physical picture, we find it useful to lift a generic four-dimensional inflating solution up to a 
five-dimensional one. This five-dimensional solution represents a pair of inflating domain wall three- 
branes and it has inhomogeneities in the additional dimension caused by the boundary potentials. On 
the other hand, initial inhomogeneities not induced by boundary sources are damped away in the linear 
regime due to the inflationary expansion and should not play any role. 

The situation is quite different in the non-linear regime, e, 3> 1, where one has to solve the full 
five-dimensional theory. To do so, we assume that the bulk scalar field has been stabilized by its 
potential and the boundary potentials allow for slow roll behaviour of the boundary scalars. Under 
these assumptions we find, in a first attempt, a simple solution by separation of variables that exhibits 
inflation. This solution represents the heterotic M-theory version of an old four-dimensional domain 
wall solution [32, |63|, recently advocated |5(| in a somewhat different approach to brane inflation. Both 
boundaries expand in a de Sitter-like manner with the Hubble parameters Hi related to the potentials 
in an unconventional way, Hi ~ |V4j|M^~ . The physical size R p h ys of the additional dimension is 
constant in time and fixed in terms of the boundary potentials. For Hi and Hi of the same order, 
H\ ~ i?2 ~ H, one has R p h ys ~ l/H. The Hubble parameter, therefore, equals the orbifold energy 
scale during inflation. We find solutions with and without an horizon at some point on the orbifold. 
The solutions without horizon require potentials with opposite sign satisfying V41 + V42 < 0. Signals 
travel from one boundary to the other in a finite time and every signal emitted somewhere in the bulk 
will eventually reach one of the boundaries. For the solutions with horizon one needs both potentials 
to be positive, V^ > 0. In this case, the two boundaries are causally decoupled. A signal emitted on 
one boundary will never reach the other one. 

In a second approach, we then find the general solution of the model without assuming separation 
of variables. This is done exploiting the similarity of our equations to those of two-dimensional dilaton 
gravity [7C]. We recover the previous inflating solution as a special case if a certain continuous parame- 
ter in the general solution is set to zero. For all other values of this parameter, however, the solution is 
non-inflating and has a collapsing orbifold. This indicates an instability of the solution which might be 
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cured by stabilizing the modulus of the additional dimension. The construction of a viable inflationary 
background in the non-linear regime is, therefore, tied to the question on how precisely such a stabi- 
lization is realized. We discuss various options and their consequences in our context. Another problem 
with the inflating solution which is made visible by its generalization is the appearance of an arbitrary 
periodic function in the solution. This function encodes the initial inhomogeneities in the additional 
dimension. Unlike in the linear case, here these inhomogeneities are not damped away. This seems to 
be in contradiction with the inflationary paradigm that all initial information should be wiped out. On 
the other hand, if sufficiently small, these inhomogeneities may lead to interesting predictions. Given 
those problems, we point out that conventional inflation in the linear regime remains a perfectly viable 
option for heterotic M-theory. For models with more than one large dimension it is the only possibility. 

Based on the results of this paper, we would like to propose three scenarios for inflation in heterotic 
M-theory. 

• The orbifold and the Calabi-Yau scales during inflation are at the specific values that at low- 
energy lead to coupling unification. In this case, the theory becomes linear for boundary potentials 
satisfying V 4i <C V^ — 6 • I0 16 GeV. At the unification point the Calabi-Yau scale and the 
fundamental Il-dimensional Planck scale are also of the order 10 16 GeV. This theory undergoes 
a transition from a pure M-theory regime at energies above 10 16 GeV (where no description in 
terms of 1 1-dimensional supergravity applies) directly to the linear regime. Inflation can then 
take place in the conventional way, presented in this paper. In this scenario, the energy density 
at the beginning of inflation is directly linked to the fundamental scales of the theory and is, in 
this sense, explained. It fits the COBE normalization V 4 ~ e 1 ' 4 6.7 • 1 16 GeV if the slow-roll 
parameter e is not too small, or, in other words, if the inflaton potential is not too flat. 

An alternative possibility is that the orbifold and Calabi-Yau scales are not be at the coupling unification 
values during inflation. Rather, they are such that the linear regime starts significantly below 10 16 GeV, 
V^ <C 10 16 GeV. In this case, there are two options. 

• Inflation takes place in the non-linear regime and is based on the corresponding solutions given in 
this paper. This option is somewhat speculative as it depends on the successful stabilization of the 
orbifold modulus. It would have very unconventional properties. These include a linear relation 
between the Hubble parameter and the potential and inhomogeneities in the orbifold direction. 
Since space-time in this context is genuinely five-dimensional, analyzing density fluctuations 
requires some care and the standard equations may not apply. 

• Non-linear inflation does not take place. This might happen for a number of reasons. For example, 
the potentials might not have the required properties, the initial conditions may not be appropriate 
or, simply, non-linear inflation might not work at all. Inflation could then start when the energy 
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density drops below Vn n and the linear regime is reached. This could be consistent with the COBE 
normalization for a very small slow-roll parameter e, that is, a very flat inflaton potential. 

2 The action 

In this section, we would like to present the five-dimensional action that we are going to use in this paper 
along with its most important properties. This includes a discussion of its origin and interpretation, 
its "vacuum" solution and the related four-dimensional effective low-energy theory that is obtained as 
a reduction on this vacuum solution. Making contact with the four-dimensional theory is particularly 
useful, in our context, whenever the relation to "conventional" four-dimensional inflation is analyzed. 
Our starting point is the five-dimensional action H 



s, = l 



2k 2 



5 UMs 

2 



R + \d a <t)d a ^ + U{(j>) 



+£ 



j=l " Jn 4 



M 



(0 



-d^d^i + u v 



(2.1) 



where the potentials are given by 



U{4>) = Ve-^ + y(</>) (2.2) 

Ui{<f>i,<f>) = TlV^ve-^ + Viifa,^ . (2.3) 

Here K5 is the five-dimensional Newton constant. Coordinates x a with indices a, j3, 7, • • • = 0, . . . , 3, 5 
are used for the five-dimensional space M$. We consider a space-time with structure M.5 = S 1 /Z% x M4 
where S 1 jZ^. is an orbifold and M.4 a smooth four-manifold. The coordinates x^ on .M4 are labelled 
by indices p,, u, p, ■ ■ ■ = 0, . . . , 3 while the remaining coordinate y = x 5 parameterizes the orbifold. It is 
chosen in the range y £ [-R, R] with the endpoints identified, where R = irp and p is the radius of the 
orbicircle. The Z2 symmetry acts as y — > —y and leaves two four-dimensional planes, characterized by 
y\ = and yi = R, fixed. These planes, separated by a distance R, are denoted by M.\ , where i = 1, 2. 



The action ( |2.1| ) describes five-dimensional gravity plus a scalar field <f> with potential U in the bulk 
coupled to four-dimensional theories on the orbifold fixed planes each carrying a scalar field <pi with 
potential Ui. The bulk fields have to be truncated in accordance with the Z2 symmetry. Specifically, 
one should require 



H-y) 


= <Kv) 


9nA-y) 


= g^u{y) 


9i&{-y) 


= -g^{y) 


955 (-y) 


= 555(2/) • 



(2.4) 



1 We have changed somewhat the notation with respect to ref. []42| to be in better accord with conventions in cosmology. 
The scalar field <f> is related to the field V of ref. |42| by V — e . The constant v was called ao there. 



Hence, 0, g^, 555 are even under the Z2 symmetry while g^ is odd. Also note that the y-derivative 
of an even field is odd and vice versa. Whereas even fields are continuous across the orbifold planes, 
odd field jump from a certain value on one side to its negative on the other side. An alternative and 
equivalent way to think about the five-dimensional space, in contrast to the orbifold picture which we 
have used so far, is the boundary picture. In this picture, the coordinate y is restricted to one half of 
the circle, y £ I = [0,R] and five-dimensional space-time is written as A4§ = I x .M4. The orbifold 
fixed planes then turn into the boundaries of this five-dimensional space. We will sometimes find it 
more convenient to work in this boundary picture. 



The potentials (|2.2| ) and (23) have been split into explicit, exponential potentials for with a height 
set by the constant v and potentials V and Vi which, at this point, are arbitrary. The — sign in ( |2.3D 
refers to the plane i = 1, the + sign to the plane i = 2. The reason for writing the potentials in this 
form will be explained shortly. It is useful to collect the mass dimensions of the various objects that we 
have introduced. In five dimensions, the Newton constant k^ has dimension three and we write 

nf = Ml . (2.5) 

Here M5 is the "fundamental" scale of the five-dimensional theory. In order to make subsequent 
equations simpler, we have pulled the Newton constant in front of the complete action. Hence, the bulk 
potentials U, V have dimension two whereas the boundary potentials Ui, Vi and the constant v have 
dimension one. 



Let us now discuss the interpretation of the action (|2.1|) -( pO| ) in terms of M-theory. More specifically, 
we will consider Hofava-Witten theory; that is, 11-dimensional supergravity on the space S 1 /Z2 x Alio 
where .M10 is a smooth ten-dimensional manifold. The two 10-dimensional orbifold fixed planes of this 
theory carry additional degrees of freedom that couple to the bulk supergravity, namely two Eg gauge 
multiplets, one on each plane. Now consider reducing this theory on a Calabi-Yau three-fold assuming 
that the radius of the Calabi-Yau space is smaller than the orbifold radius. For the present values 
of these radii, such a relation is suggested by coupling constant unification pq| . We then arrive at a 
sensible five-dimensional theory on the space-time M5 = S /Z2 x .M4, where the two four-dimensional 
fixed planes of the orbifold result from the original 10-dimensional planes. For the standard embedding 
of the spin connection into one of the E% gauge groups, this effective action has been computed in 



ref. [42, 52, |43|| . The generalization to non-standard embedding has been described in ref. |5j]. It 
turns out that the bulk theory is a five-dimensional N = 1 gauged supergravity coupled to vector- 
and hypermultiplets. This bulk theory is coupled to two four-dimensional jV = 1 theories that reside 
on the now four-dimensional orbifold planes. More specifically, these boundary theories contain gauge 
multiplets as well as chiral (gauge matter) multiplets. Upon appropriate reduction on the orbifold to 
four dimensions (in a way to be specified below), one obtains a candidate for a "realistic" jV = 1 
supergravity theory with the observable sector coming from one plane and the hidden sector from the 
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other. 

The action ( |2.lD is a "universal" version of this five-dimensional effective theory truncated to the 
field content that seems essential to discuss cosmology. Let us try to make this statement more precise 
by explaining the meaning of the various objects in fl2.1| )-(pli|) in terms of the underlying 11-dimensional 
theory. In terms of the 11-dimensional Newton constant k and the Calabi-Yau volume vcy the five- 
dimensional Newton constant is given by 

4 = — ■ (2-6) 

The bulk field 4> is simply the modulus associated with the Calabi-Yau volume such that the "physical" 
volume is given by e^t>cY- Clearly, the Calabi-Yau space has, in general, many more moduli associated 
with its shape and complex structure. For simplicity, we have kept the volume modulus only, since it 
is the geometrical modulus common to all Calabi-Yau spaces. It is in this sense that we are referring 
to the action as universal. Also, for simplicity, we have dropped a number of other bulk fields such 
as the axions associated with <p an d the vector fields in the vector multiplets. All these fields can be 
consistently set to zero in the full equations of motion. Therefore, our solutions will be solutions of the 
complete action of heterotic M-theory as well. We have, however, kept a feature of the action that arises 
from the gauging of the bulk supergravity and is essential to discuss cosmology; namely, the explicit 
bulk potential for in eq. (|2.2|). The constant v in this potential is given by 



TV / k \ 2 / 3 n . 

" = -7s(*f) j* (2J) 

where n is an instanton number related to the tangent space and gauge instantons in the internal 
Calabi-Yau space. We have added another potential V{<p) in the bulk which one expects to arise from 
non-perturbative effects like internally wrapped membranes. The explicit computation of this potential 
for M-theory is not very well understood at present (see however [J55], |5(|). It is, however, clear that 
such a potential is eventually needed to stabilize the moduli. Given the lack of theoretical knowledge, 
we will not assume any specific origin for V but, rather, try to specify what properties it should be 
required to have from a cosmological viewpoint. 

Let us now move on to the four-dimensional theories on the orbifold planes and their M-theory 
origin. As already mentioned, these planes arise directly from the 10-dimensional planes of Hofava- 
Witten theory upon reduction on the Calabi-Yau space. On each of these planes, we have introduced 
a scalar field <pi with potential Ui that represents the scalar partners of the matter fields. In fact, for 
non-standard embedding such matter fields generally arise on both planes. Of course, restricting to 
one scalar field on each plane is a tremendous simplification from the particle physics point of view. 
Cosmologically, however, it seems reasonable to choose such a model with one "candidate inflaton" on 
each orbifold plane, especially in a first study of inflation in such models. Although the potentials Vi 
in eq. ( p. 3D are known in principle for a given Calabi-Yau compactification, here we will not attempt 
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to be more specific about their form. Instead of going into such detailed questions of model building, 
we will assume that they have the cosmologically desired properties. There is another, explicit part 
of the boundary potential in eq. ( |2.3| ) which depends on the projection of the bulk field eft onto the 
orbifold planes. These potentials originate directly from Hofava-Witten theory and are related to the 
explicit bulk potential for 4>. Note that, in particular, the height is set by the same constant v, defined 
in eq. ( |2.7D . As we will see, they support a three-brane domain wall solution of the five-dimensional 
theory. In the following, we, therefore, refer to them as domain wall potentials. Of course, there will be 
gauge multiplets on the orbifold planes as well which we have not written in ( |2.1| ). They could play a 
role in cosmology via gaugino condensation. However, we will not consider this explicitly in the present 
paper. To summarize, the action ( |2.1| )-( |2~3| ) can be viewed as part of the five-dimensional effective 
theory of heterotic M-theory and it contains the basic cosmologically relevant ingredients of this theory. 



The M-theory context is not obligatory here. Instead, the action (^Tj) could be viewed in the 
general context of theories with large additional dimensions where the standard model arises from a 
brane worldvolume. In fact, in this context, our action is about the simplest appropriate to study 
inflationary cosmology. We have one "large" dimension (the orbifold) and two four-dimensional brane- 
like theories corresponding to an observable and a hidden sector. In addition to bulk gravity, we have a 
bulk scalar field with potential and a scalar field with potential on each plane; in other words, a minimal 
setting for inflation. The explicit potentials in (|2,2j), ( |2.3| ) that originate from M-theory can always be 
switched off by setting v = 0, if desired. Our results are, therefore, not limited by their presence. What 
does limit our results is the presence of only one additional dimension, something that is appropriate for 
heterotic M-theory but not necessarily otherwise. We will comment on the modification of our results 
in the case of two or more large dimensions as we proceed. 



For later reference, let us collect the equations of motion derived from the action ( |2.lD . We have the 
Einstein equation 



R a p - -^g a pR = T af3 + ^2 T al s (y ~ Vi) ( 2 - 

8=1 



with the bulk and boundary energy momentum tensors 



T a p = -\ (daWtf - ^g a p(d0)A + X -g a(i U (2.9) 

r$ = -^g^ 2 (d^d^i - ^u(d^i) 2 - g^U?j (2.10) 

r$ = o, rf = o. (2.ii) 



For the scalar fields we find 

2 
□ _,A _ Fl.TT _ ^~V 2 



9^-555 1/2 ^^^(y- yi ) = (2.12) 



i=i 



□4^-^V- = (2.13) 



where Ds and D 4 are the Laplacians associated with the five-dimensional metric <7 Q/ g and its four- 
dimensional part g^y projected onto one of the orbifold planes. 

Which types of solutions to these equations are we interested in? For a cosmological solution, 
one would like to have a three-dimensional maximally symmetric subspace which, for simplicity, we 
take to be flat. Studying open and closed universes would be clearly interesting as well. Later, we 
will distinguish between a linear and a non-linear case. While the generalization to include spatial 
curvature is straightforward in the linear case, the non-linear case is significantly more complicated. 
Still, we do not expect our main conclusions to depend on the choice of spatial curvature. Clearly, the 
maximally symmetric subspace cannot contain the orbifold. Hence, the solutions are independent of 
x = (x 1 , x 2 , x 3 ). They can, however, depend on the time r = x° as well as on the orbifold coordinate y. 
As we will see in a moment, the presence of fields and potentials on the boundaries in fact forces us to 
consider orbifold dependence. Therefore, we start with the ansatz 

(2.14) 
(2.15) 
(2.16) 
(2.17) 
(2.18) 
(2.19) 

Here a is the scale factor of the three-dimensional universe and (3 is an additional scale factor that 
measures the orbifold size. In some cases, we will choose the conformal gauge v = (3 in the above 
metric. Note that we can choose this gauge in the (r, y) subspace in such a way that the boundaries 
are mapped to hypersurfaces y = const. Usually, then, we can shift the boundaries back to y = 0, R so 
that our conventions for the coordinate system remain intact. A special case occurs if a boundary has 
been "mapped to infinity" due to a singularity in the reparameterization. To be able to deal with this 
case, we keep v arbitrary in the subsequent formulae. 

The equations of motion for the above ansatz are given by 



dsl 


= 


-e 2v dr 2 + e 2a d^ + e^dy 2 


V 


= 


v(t, y) 


a 


= 


a(r,y) 


(3 


= 


P(r,y) 


4> 


= 


4>{r,y) 


(j).; 


= 


Mr) ■ 



2,e- 2u {a 2 + a/3) - 3e" 2/ V - ol + 2a' 



2 

V V + V 2 V 2 + -U + -e~? Y 
4 Y 4 r 2 2 *-" 
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\e- 2v $ + U t 



5{y- yi ) (2.20) 



e~ Zv (2a + + (3d + 2/3 - 2z>)d + {(3 - u)/3) 



e -W{2a" + v" + (3a' - 2(3' + 2v')o! + (i/ - £>' 



r^ 2+ k 2V + ^-r 



-/3 



E 



1 -2v \2 TT 

2 e & - V i 



6(y-Vi) (2.21) 



3e- 2i/ (d - va + 2d 2 ) - 3e" 2/3 (a /2 + i/V) 
3(d' + aa' — av' — 0a') 

e- 2u (4> + (3d + $ - v)4>) - e- 2l3 {4>" + (3a' - + i/)<//) 

4> t + (3d-i>)^ i + e 2u d <j> y t 



1 



1 



1. 



-e 2 > 2 -^e- 2 V 2 + -£/ (2.22) 



o 



(2.23) 

d^U - e~P Y, d <P U Av ~ Vi) ( 2 - 24 ) 

(2.25) 



i=i 



where the dot and the prime denote the derivatives with respect to r and y. Alternatively, one can 
formulate these equations in the boundary picture. Then all 5 function terms in the above equations 
have to be dropped and are replaced by the following boundary conditions 



(2.26) 
(2.27) 
(2.28) 



e"V 1 


1 
= ^12 


^ -2v 12 i TT 

2 e 9% + Ui 


y 


^v U, 


1 

= T 12 


--e 0i + ?7i 


--"*' U, 


2 L 


^4=^ • 





y=yi 



Here the upper (lower) sign applies to the first boundary at y = y\ (the second boundary at y = 2/2)- The 
(5-function terms, or equivalently the right hand sides of the boundary conditions, are non-vanishing 
if there is any kinetic or potential energy on the boundaries. Hence, in this most interesting case, the 
solution is necessarily inhomogeneous in the orbifold coordinate y. 

For the application to cosmology, one is also interested in the four-dimensional effective action 
of fl2.1|) that is, roughly speaking, valid when all energy scales are smaller than the orbifold scale 1/R. 
Normally, this action could be derived by a simple truncation where all the bulk fields are taken to be 
independent of the orbifold coordinate. While, in our case, this gives the correct answer to lowest order, 
it neglects higher order corrections that can be relevant. These corrections appear because, strictly, we 
are not allowed to take the bulk fields independent of the orbifold coordinate, thereby neglecting all 
contributions from Kaluza-Klein modes. Instead, as we have just seen, every non-vanishing term in the 
boundary actions produces an orbifold dependence that needs to be integrated out and, typically, leads 
to corrections to the effective action. Let us explain how this works in our case, taking into account 



corrections up to the first non-trivial (linear) order. For a more detailed account see [57]. First, we 
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split the bulk fields as 

9a/3 = 9a/3 + 9a/3 (2-29) 

<j> = + (i> (2.30) 

into their orbifold average plus an orbif old-dependent variation. Specifically, we have defined 

9af3 =< 9aj3 >5 , <t> =< 4> >5 (2-31) 

where < . >5 denotes the average in the orbifold direction. Hence, the average of the variations vanishes; 
that is, 

< ~9 a p >5= o , < 4> >5= . (2.32) 

The averaged fields are in one to one correspondence with the low energy moduli fields. Aside from the 
modulus 4>, we have the orbifold modulus (3 and the four-dimensional Einstein-frame metric g^v The 
latter two quantities are related to the averaged metric by 

e W = 555 (2.33) 

9A„u = e^g^. (2.34) 

Note that we have no graviphoton zero mode since g^ is odd under Z 2 and, hence, g^ = 0. After 
expanding to linear order in the variations, the equations of motion can be decomposed into an averaged 
part and a part for the variations. The latter equations have the following solutions in the boundary 
picture R 



9 a p = e 



2 ?R 



P 1 (z)S^+P 2 (z)S^ (2.:W) 



\e?R [PxWd+Ux + P2(z)d <t> U 2 ] (2.36) 



where 



S% = T§ ~ \ 9a ^ 5 T§ (2.37) 

is the modified boundary stress energy. The bar indicates that the corresponding expressions are under- 
stood with the bulk fields replaced by their zero modes. Furthermore, we have defined the polynomials [j 

P 1 (z) = — z 2 + z--, P 2 (z) = — z 2 + - (2.38) 

v ' 2 3 w 2 6 v ; 



2 This solution applies to the part of the boundary stress energy which is homogeneous and constant in time on scales of 
the order R. Slowly varying potential energy, the most important case for this paper, typically satisfies this requirement. 
For boundary processes on scales smaller than R the associated bulk fields are suppressed and decay exponentially away 
from the boundary |5q ]. 

3 These polynomials are characterized by the properties < Pi(z) >5 = < ft(^) >5= 0, P{(0) — 1, -Pi'(l) = 0, ^(O) = 
andP^l) = -1. 
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which depend on the normalized orbifold coordinate 



R 



G[0,1]. 



(2.39) 



This solution shows that each energy source on the boundary leads to a certain, generally quadratic, 
z-dependent variation of the bulk fields across the orbifold. This corresponds to coherent excitations 
of the Kaluzu-Klein modes in the orbifold direction. Let us consider an explicit example. We would 
like to determine the z-dependence of the bulk fields that arises from the explicit boundary potentials 
in ( |2.3D proportional to v. Inserting the corresponding part of the boundary stress energy into eq. fl2.35|) 
and ( 2.36| ) we find 



dSr 



(g a /3 + g af 3)dx a dx 13 

( 1 



1 + ^0 e- p g 4 ^dx^dx v + 1 + 



4~ 



e 2 ?dy 2 



-2e?-+< 



where 



£DW 



vR 

71 



(2.40) 
(2.41) 
(2.42) 

(2.43) 



Note that the explicit potentials have the same height but opposite sign on the two boundaries. There- 
fore, the z-dependent part of the solution is proportional to the difference (Pi —P^jiz) of the polynomi- 
als (|2.3§| ) and, hence, is linear in z. This is to be contrasted to the general case of unrelated boundary 
potentials which leads to a quadratic variation. For constant moduli g^v, J3 and (f>, eq. ( 2.40| )— ( 2.43 ) 
represent the linearized version of an exact BPS domain wall (three-brane) solution of heterotic M- 
theory [] in five dimensions that was found in ref. [42]. More precisely, the solution represents a pair 
of domain walls stretched across the orbifold planes and it can be viewed as the "vacuum" solution of 
the theory. At the same time, it is the five-dimensional version of Witten's 11-dimensional linearized 
background [pS(| . In this way, the four-dimensional orbifold planes are identified with three-brane world- 
volumes which carry the observable (and hidden) low energy fields as their zero modes. In this sense, 
our picture offers a concrete realization of the general idea that the world arises from the worldvolume 
of a brane. Observe that the size of the correction (2.40) is set by the dimensionless quantity £dw 



in eq. ( 2.43 ), which is just the product of the boundary potential (measured in units of k 5 = M^) 



times the size of the orbifold dimension. Therefore, the linearized approximation that led us to the 
solution (2.35), ( 2.36J ) is only sensible as long as |cdw| ^ 1- We will discuss this in more detail in the 
next section. 



4 The exact solution of ref. |42[ has the form ds% = a^Hri^dx* 1 dx v + b^H^dy 2 



bgH 3 where H = cq — |eow(z — o) 



and ao, bo and Co are constants. In fact, it also constitutes an exact solution of the action (2.1). Upon linearizing in edw, 
setting g4/iv ~ f]p,v in (2.40) and appropriately matching the moduli this coincides with (2.4C )— (J2.42|) . 
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By inserting the solution ( 2.29|) -( p.36|) into the action (^]|) and promoting the moduli g^ v , J3 and 
4> to four-dimensional fields, we can now compute the first order corrected action of the zero modes. 
We introduce the fields 



S = e^ , 



(2.44) 



which are just the real parts of the ordinary S— and T-moduli []. Note that S measures the volume of 
the internal Calabi-Yau space in units of vcy whereas T measures the size of the orbifold in units of R. 
We also introduce the four-dimensional Newton constant 



Gn 



1 



16vri? 16irRMl 



and the rescaled boundary fields 



a 



mi 



Then the four-dimensional effective action is given by 



Si 



1 



16ttGn 



\/ 3 34 



M 4 



1 



Ra + -T- 2 d a Td»T + -S~ 2 d u Sd^S 



r i 2 

J M 4 l ~7 



with the "Kahler metrics" 



i=\ 



K, 



3 , epw 
IT 25 



and the four-dimensional potential 

2 



* = 2^> + H £ (4 V « ± d * V «) + T- RM * V 

i=X i=l ^ ' 



where the boundary potentials V^ normalized to four mass dimensions are defined by 



V Ai = MgVi . 



(2.45) 



(2.46) 



(2.47) 



(2.48) 



(2.49) 



(2.50) 



Note that the Kahler metrics and the potential have corrections linear in edw that originate from 
the domain wall. Had we performed a simple truncation of the five-dimensional theory by taking the 
bulk fields independent of the orbifold we would have missed these corrections. What happened to the 
explicit potentials proportional to v in the five-dimensional action, eq. ( p.2[ ) and ([2.3])? These potentials 
were actually responsible for the existence of the domain wall solution in the first place. Performing 



5 The imaginary parts are absent because we have omitted the corresponding fields in the five-dimensional action ( |2.]| ) 
for simplicity. The complete reduction can be found in ref. [43] ■ 
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a reduction on this solution, they are canceled so that there is no remnant potential for the S and T 



moduli in four dimensions. Given the form of our five-dimensional action (2.1)- (pO| ), the bulk potential 
V only depends on S but not on T. The four-dimensional effective potential fl2.49| ) shows that, in this 
case, the modulus T cannot be stabilized. For the solutions based on the four-dimensional effective 
action we, therefore, have to slightly generalize our setup and assume that V is a function of the orbifold 
modulus T as well. 

An alternative way to deal with the orbifold dependence of bulk fields in a reduction to four di- 
mensions is to keep all the Kaluza-Klein modes instead of integrating them out. Recently, this has 



systematically been carried out in ref. [59|. In such a four-dimensional action, the zero mode part does 
not have the above corrections. Instead, there exist terms in the action which couple Kaluza-Klein 
modes linearly to zero mode fields. Hence the Kaluza-Klein modes cannot be set to zero consistently. 
In fact, solving for these Kaluza-Klein modes would lead to the Fourier decomposition of our orbifold- 
dependent corrections fl2.35| ), ( 2.36 ). 



Besides being useful to calculate the four-dimensional action, the method described above can also 
be used to find approximate solutions of the five-dimensional theory. This is actually the main reason 
why we have presented it here in some detail. Concretely, suppose one has found a solution of the 
four-dimensional theory with action (|2.47l ). Then, by inserting this solution into eq. ( p.29|) - (|2.36| ), we 
can simply "lift it up" to obtain a solution of the five-dimensional theory. Of course, for this solution to 
be sensible, the linearized approximation that led us to eq. (2.35|), (|2.36) should be valid. The condition 



for that to be true will be discussed in the next section in some detail. 



3 Types of inflation and scales 



Which type of inflation for the action (2.1) do we want to consider? First of all, in this paper we are 



interested in "conventional" potential-driven inflation rather than in a pre-big-bang-type scenario [6C|. 



Some solutions of five-dimensional heterotic M-theory that might provide a basis for a pre-big-bang 



scenario have been found in ref. pi| , 32, p3|| . Given that we want to focus on potential energy, there are 
two obvious options. 

• bulk potential energy: The potential energy is provided by the bulk potential V and the field (p is 
the inflaton. This can also be called "modular inflation". 

• boundary potential energy: The potential energy is due to the boundary potentials V% (one of 
them or both) and the fields 4n are the inflatons. This can be called "matter field inflation". 

Clearly, in general, one could also have a mixture of both types. In this paper, we concentrate on 
the second option of boundary potential energy, which appears to be more interesting for a number of 
reasons. Most importantly, the presence of the boundaries is the truely new ingredient in the action ( |2.1| ) . 
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Inflation from the boundary also seems to be in accord with the current mainstream in four-dimensional 
inflationary model-building |pl| . On the other hand, modular inflation faces a number of problems 
associated with the steepness of typical non-perturbative moduli potentials p4| , |65| ], From the viewpoint 
of the recently proposed models [^(J with a very low "fundamental" scale M5 ~ TeV, bulk inflation 
might not be desirable because reheating from a bulk field might leave the additional dimension too 



inhomogeneous to be consistent with standard cosmology [47]. None of these arguments, of course, 



disproves bulk inflation and it might still be an interesting option. This will be discussed elsewhere. 

Let us, from now on, concentrate on the case of potential energy on the boundary. Our assumption 
is, hence, that the bulk potential V{<p) is zero or negligible for all field values 4> relevant during inflation 
(which does not mean, however, that V is identical to zero). The potential energy is then dominated 
by the boundary potentials Vi and the four-dimensional effective potential is given by 

1 2 
Vi ~ ^2 E V " + °( e Dw) , V Ai = M%Vi . (3.1) 

In the previous section, we have given a method to obtain approximate five-dimensional solutions 
from solutions of the four-dimensional effective action. Therefore, if we can find a four-dimensional 
inflationary solution (which we can if V4 satisfies the usual slow roll conditions), we can lift it up to a 
solution of the D = 5 action. When is this effective four-dimensional approach sensible? In the domain 
wall example of the previous section we have seen that the condition 

\v\R 
|cdw| = M=- < 1 (3.2) 

which involves the height v of the domain wall potentials should be satisfied in order for the linearized 
solution to be valid. In the context of heterotic M-theory, the relation ( |3.2[ ) is required anyway since 
the formulation of the theory is not well known beyond the linear level in enw- in an y case, (|3.2[) should 
be satisfied for the effective four-dimensional action fl2.47| ) to be sensible. 

Inspection of the general linearized solution ( 2.35J ), ( |2.36| ) shows that we should have a similar 
condition for the boundary potentials Vi] namely, for 

e * = ViR=Y nf (3 - 3) 

we should have 

M « 1. (3.4) 

What is the interpretation of these conditions? We have generally seen that boundary potentials (and 
any other form of boundary energy) lead to inhomogeneities in the additional dimension or, in other 
words, excite Kaluza-Klein modes associated with this dimension. If |e$| <C 1 this can be described 
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in a linearized approximation, either by integrating out the Kaluza-Klein modes at the linear level as 
we have done to arrive at our four-dimensional action for the zero modes, or by keeping the linearized 
Kaluza-Klein modes in the four-dimensional action as in ref. EiJ. If, on the other hand, |ej| 3> 1 this 
linearized approximation breaks down. Then the simplest possibility is probably to work with the full 
five-dimensional action. Hence, we distinguish two cases. 

• |ej| <C 1, the linear regime: We can solve the four-dimensional effective action ( |2.47| ) for the zero 
modes and lift up the solutions to five dimensions using the results of the previous section. We 
will discuss this case in the following section. 



• |ej| 3> 1, the non-linear regime: We should solve the full five-dimensional theory fl2.1|) . This case 
will be considered in section 5 and 6. 

What seems surprising about this criterion, at first, is that the quantities e« are linear in the size 
R of the additional dimension. The larger the additional dimension, the earlier the system enters the 
non-linear regime. This can be understood as follows. Consider the one-dimensional space associated 
with the additional dimension. The orbifold planes appear as points sources in this space. On the linear 
level, the fields generated by these point sources are roughly described by the one-dimensional Green's 
function which is simply the linear function \y\. The proportionality of 6j to the orbifold size R simply 
reflects this linear increase in y of the one-dimensional Green's function. This picture also suggests how 
the above criterion should be modified for more than one additional dimension. For example, for two 
additional dimensions the Green's function is a logarithm and, hence, roughly a constant. Therefore, €j 
should be independent of the size of the additional dimensions. Generally, for d additional dimensions, 
one expects the corresponding parameter e^ (suppressing the index i) to be given by 

€d ~ MlXJR*-* { ' 

where M^ + d is the fundamental scale of the 4 + d-dimensional theory. Since one would generally require 
that V4 <C Mf +d and R^ 1 <C M 4+ ^ to have a field-theoretical description, it follows that |e<j| <C 1 as 
long as d > 1. We see, therefore, that the case of one additional dimension, which is the relevant one 
for heterotic M-theory, is special in that it is the only case where inflation might take place in the 
non-linear regime. 

We would like to be somewhat more explicit about scales in order to get a feeling for when the 
above linearity criterion might be satisfied. In the Hofava-Witten context, one reference point is the 



"physical point" at which the gauge and gravitational couplings are matched [28]. At this point, one 
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finds to lowest order 



K 



~ 2/9 ~ x4-10 16 GeV 



v~y /6 ~ x 2 • 10 16 GeV (3.6) 

R- 1 ~ x 3 2 • 10 15 GeV 

for the energy scales associated with the 1 1-dimensional theory, the Calabi-Yau space and the orbifold. 
Here x is a quantity of order one which depends on the shape of the Calabi-Yau space and parameterizes 
our ignorance of the precise relation between the Calabi-Yau volume and the unification scale. Using 
eqs. (|2.5|) and (|2.6|), we find for the mass scale of the five-dimensional theory 



M 5 ~x2-10 17 GeV. (3.7) 

The linearity criterion (|3.4|) then translates into 

Vlt « x 3/2 6 • 10 16 GeV . (3.8) 

Hence, if during inflation the Calabi-Yau volume and the orbifold radius are at the values that lead 
to coupling constant unification, the Kaluza-Klein modes behave linearly as long as the boundary 
potentials satisfy the bound ( |3.8| ). This has to be compared with the COBE normalization which 
implies for the (four-dimensional) inflationary potential that 

y 4 1/4 ~e 1/4 6.7-10 16 GeV (3.9) 

where 

- ^ (3.10) 



2G N Vi 

is the usual slow roll parameter. The prime denotes the derivative with respect to the inflaton. Compar- 
ison of eq. ( |3.9D and eq. (|3.8| ) shows an interesting coincidence of scales suggesting that inflation might 
start just when the theory leaves the non-linear regime. The bound ( |3.8j ) is, however, very close to the 
physical orbifold and Calabi-Yau scales in ([O]) . Beyond energies of v CY the effective five-dimensional 
theory breaks down. Also, beyond energies of k~ 2 ' 9 the description via 1 1-dimensional supergravity 
is no longer viable. Clearly, beyond those scales our analysis does not apply and there is no sense to 
talk about a non-linear regime as defined above. Most likely, therefore, the theory undergoes a tran- 
sition from an M-theory regime directly into the linear regime. Of course, the "coincidence" of scales 
is due to the fact that the fundamental scales of our theory are much closer to the scale suggested by 
COBE than, for example, the Planck scale. To make this observation really meaningful one has to more 
closely analyze the transition to the linear regime, something that is beyond the scope of this paper. 
In any case, all these statements relate to the "physical point" associated with the present values of 
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the Calabi-Yau volume and the orbifold radius. Clearly, those values could have been different in the 
early universe, so that we really do not know in which regime inflation took place. Consequently, in the 
context of Hofava-Witten theory, as well as in a wider context, we should investigate both the linear 
and the non-linear possibility. 

Finally, we would like to discuss which part of the space should be allowed to inflate from a "phe- 
nomenological" point of view. There are two options. First, while the usual three-dimensional space 
inflates the additional dimension is basically fixed. Secondly, both the three-dimensional space and the 
additional dimension inflate. While there is obviously no problem with the first option, inflating the full 
space might lead to an additional dimension that is too large. In particular, inflating it for the full 55 or 
so e-folds leaves one with a radius R that would typically be some 24 orders of magnitude larger than 
the fundamental length scale M§ of the theory. This is clearly unacceptable. It might be acceptable 
to inflate the additional dimension for a short period (depending on the precise value of the scales) 
and then stabilize it to its low energy value. This stabilization of an inflating modulus might, however, 
be hard to achieve theoretically. In this paper, we therefore favor the first option of a non-inflating 
additional dimension. 

4 Linear case: The inflating domain wall 

In this section, we would like to construct inflationary solutions in the linear case 

M « 1 • (4.1) 

Following the general method presented in section 2, we should first find an inflating solution of the 
four-dimensional effective theory ( 2.47| ). The existence of such a solution depends, of course, crucially 



on the properties of this theory and its potential in particular. Rather than going into model building, 
we will simply assume suitable properties. Since we are interested in boundary inflation, we assume 



that the moduli S and T have been stabilized by the four-dimensional effective bulk potential (2.49). 
As discussed earlier, this requires a non-perturbative bulk potential V that depends on S and T. To 
simplify formulae, we choose the coordinate Calabi-Yau volume vqy and the coordinate orbifold radius 
R such that this stable point is at 

= lnS = O, /3 = lnT = 0. (4.2) 

Furthermore, since we would like the vacuum energy to be dominated by the boundary potentials, we 
assume that the bulk potential vanishes at this point; that is 

V{<j> = 0,j3 = 0) = 0. (4.3) 
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Finally, we need our candidate inflatons, the boundary fields Cj, to be slowly rolling. This requires the 
inequalities 



(^ 4 ) 2 <<l5 Oc^CjVa c 1 



G N V 2 



GmVa 



NV4 



(4.4) 



to be satisfied for all i,j = 1, 2. Then, starting with the usual four-dimensional metric 

ds\ = g^dx^dx" = -dr 2 + e 2a4 dx 2 (4.5) 

for a spatially flat Robertson- Walker universe with scale factor 04, the four-dimensional equations of 
motion from the action (2.47]) reduce to 



2 _ .2 87rgjy 
-H = a 4 = — - — 14 



C, 



1 



9c^4 • 



The usual inflating solution is then 



a 4 = H(t - t ) 



(4.6) 
(4.7) 

(4.8) 



with an arbitrary integration constant To. The slow roll motion Ci{T) of the inflatons can be obtained 
from eq. ( [4.7| ) once an explicit potential has been specified. Note that with the conventional relation fl4.6|) 
between the potential and the Hubble parameter and using eq. (2.45), the quantities e« can be written 



as €i ~ H 2 R 2 . The condition £j C 1 then corresponds to the "intuitive" criterion that the Hubble 



parameter H should be smaller than the mass R~ l of the first Kaluza-Klein excitation. In ref. [49] 
this criterion has been used to constrain the Hubble parameter during inflation for TeV-scale gravity 
models. 

We can now lift this solution up to five dimensions using the formulae of section 2, in particular 
eq. (|2~29|)-(]2~36D. We find | 

-e 2 "dT 2 + e 2a d^ + e 20 dy 2 (4.9) 

H (r - to) + a (4.10) 

a (4.11) 

4a (4.12) 

^^ 4> (4.13) 

D A subtle point has to be taken into account if one wants to explicitly verify that this solution satisfies the five- 



dsj 
a 
v 
P 



dimensional equations of motion ( 2.20 )-( 2.25) to linear order. For the underlying four-dimensional solution, we have 

The effective four-dimensional 
1) = and V(T = 1) = 0, 



explicitly assumed that the orbifold modulus T — e p has been stabilized at T 



potential (2.49) shows that this requires aT dependence of the bulk potential V . Since OtV^T 

^2 



one concludes from eq. ( 2.4S ) that dpV(T — 1) = R 1 2»=i ^ + O(cdw)- Moreover, in the five-dimensional equations of 
motion we have not considered a /3 dependence of the bulk potential V . To incorporate such a case, the potential U in 



cq. ( 2.22) has to be replaced by U + dp V . Using this modification and the above expression for dp V one can indeed verify 
that the five-dimensional equations of motions are satisfied. 



19 



where the orbifold dependent corrections a and 4> are given by 
a 



( z ~l)- ^sEfl(^(^ = O.Cf 4 (r)) (4.14) 



2 



= -2e^(z-^\+^Y j P i {z)d 4> V A S = ^C i {T)). (4.15) 

\ / 5 j = i 

Recall that P%{z) are quadratic polynomials in the normalized orbifold coordinate z = y/R £ [0, 1] that 



have been defined in eq. ( 2.38 ). In the metric ( ff-9 ) we should, of course, only consider terms linear in 
a and 0, in accordance with our approximation. 

Let us discuss the form of this solution. First we note that, neglecting the very mild time dependence 
introduced by the slow roll of the inflatons for the moment, the solution separates into a time-dependent 
and an orbifold-dependent part. The time-dependent part just corresponds to the inflationary expan- 
sion of the three-dimensional universe. This expansion does not occur just on one boundary, as one 
might naively expect, but uniformly across the whole orbifold. To discuss the corrections, let us con- 
centrate on the scale factor a. Although it expands at the same Hubble rate H everywhere across 



the orbifold, its actual value depends on the orbifold point as specified in eq. ( 4.14 ). The first term 
in this equation is the familiar linear contribution from the domain wall proportional to €dw- The 
second term arises from the boundary potentials and is proportional to ej, as expected. It has a mild 
time-dependence through the slow-roll change of the potentials. 

The bottom line of this section is, that the problem of finding inflationary backgrounds in the linear 
regime can be adequately approached in the four-dimensional effective action obtained by integrating 
out the Kaluza-Klein modes. For our simple model, this action is given in ( |2.47 ). More realistic four- 



dimensional effective actions from Hofava-Witten theory can be found in p^, 58, B8|, 54, 39 1. The 
full five-dimensional solution is then obtained by lifting the four-dimensional solution up, using the 



correspondence established in section 2. This leads to the corrections (4.14) and (|4.15|) corresponding 
to Kaluza-Klein modes that are coherently excited by the non-vanishing sources on the orbifold planes. 
One might also worry about other excitations of the Kaluza-Klein modes unrelated to the orbifold 
sources, such as remnants from an initial state. This could be described by adding the tower of Kaluza- 
Klein modes to the four-dimensional effective action ( [2. 47 ). Since we have integrated out the orbifold 
sources, those modes would be free source-less particles with masses n/R, where n is an integer. During 
inflation, oscillations of these modes are simply damped away by the expansion. The condition for this 
to happen efficiently coincides with our linearity criterion |ej| <C 1 and is, hence, satisfied. In the 
linear regime, the only relevant excitations of Kaluza-Klein modes after a short period of inflation are, 
therefore, the ones caused by the orbifold sources computed above. As we will see, this changes in the 
non-linear regime where |e^| 3> 1. 

We would briefly like to mention some generalizations. It is clear that the above method can be 
applied to other types of four-dimensional cosmological solutions, for example to a preheating solu- 
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tion with the energy density dominated by coherent oscillations or to a radiation dominated solution, 
straightforwardly. Basically, all one has to do is to replace the potentials in eq. ( [4.14 ) by the appropriate 



energy density. What about the case of more than one large dimension? In such a case, the variation 
across the orbifold would not increase as a polynomial any more, as we have seen earlier. Instead, for 
two additional dimensions we expect a logarithmic behaviour and for more than two dimensions a power 
law fall-off. Common to all these cases is, however, that energy density on the four-dimensional plane 
coherently excites the bulk modes. 

5 Non— linear case: A solution by separation of variables 

We would now like to study the non-linear case; that is, we assume 

M » 1 • (5.1) 



In this case, we should solve the full five-dimensional equations of motion given in eq. ( 2.20[ )-( ^.28[ ). We 



will use the boundary picture to do this. In general, we have two types of potentials on the boundaries, 
namely the domain wall potentials and the potentials Vi, corresponding to the two terms in ( |2.3j ). We 
have already stated that eDW> the dimensionless quantity that measures the strength of the domain 
wall corrections, should be smaller than one in order for the M-theory description via supergravity 



to be valid. The condition (5.1), therefore, states that the potentials Vi will be dominating over the 
domain wall. To simplify our problem, we will therefore neglect the domain wall potentials. Certainly, 
there will be an intermediate region between the non-linear and linear regime where both potentials 
are significant. It will, however, be very difficult to find explicit solutions in this regime. We, therefore, 



concentrate on the case (5.1). As a further simplification, let us assume that the Calabi-Yau volume 



modulus <j) has been stabilized by the bulk potential V; that is 

4> = const . (5-2) 

Of course, we have to be careful that this assumption is consistent with the boundary condition on (f>, 
eq. ( |2.28| ) . We have already neglected the first term in this condition which originates from the domain 
wall. The second term is related to the boundary potentials Vi and vanishes if those are taken to be 
independent of <p. We will assume this in the following. In accordance with our general assumption of 
boundary inflation, the potential energy from the bulk potential should be negligible, 

V(<f>) ~ . (5.3) 

Finally, we assume that the boundary potentials Vi(4>i) are suitable slow-roll potentials so that the 
boundary fields (pi act as the inflatons. Practically, this means that we treat Vi simply as constants. 
We recall from section 2 that the metric has the form 

dsl = -e 2 »dr 2 + e 2a dx 2 + e 2f3 dy 2 (5.4) 
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where a and /3 are the scale factors of the three-dimensional universe and the orbifold respectively. We 
also remark that, from eq. ( p.26| ), ( |2.27j ), the boundary conditions take the form 



Even though the equations of motion ( 2.20[) -( p,28 ) have now considerably simplified, they are still 



not easily soluble. The reason is, of course, that we are dealing with partial differential equations as 
opposed to ordinary ones that one usually encounters in cosmology. The simplest strategy to solve 
partial differential equations is separation of variables and this is what we are going to do next. The 
general solution will be given in the following section. 

To simplify the problem we first choose the conformal gauge v = (3. We will assume for the time 
being that the coordinate transformation that led to conformal gauge leaves the boundaries at finite 
values of the coordinate y. In this case, we can restore the conventions for our coordinate system by 
shifting the boundaries back to y = 0, R. We are, then, looking for all solutions of fl2.20 )-( 2.2S| ) (subject 
to the above assumptions) consistent with the separation ansatz 

a = a (r) + a 5 (y) (5.6) 

/3 = A)(r) +&(</)• (5-7) 

The general solution to the equations of motion ( p.20[) -( p.23[ ) in the boundary picture is 

a = K(y±r)+A (5.8) 

j3 = K((l-K)y±kr)+B (5.9) 



where K, K, A and B are integration constants. We still have to apply the boundary conditions ( |5.5[) . 
This leads to 

*->(-$)• *-■ —(-") 

while A remains arbitrary. For the arguments of the logarithms to be positive, we have to further 
demand that V\ and Vi have opposite signs such that V\ + V% < 0. It is not yet clear, whether these 
restrictions on Vi are general or whether they are related to our choice of the coordinate system. We 
have assumed that in conformal gauge the boundaries are at finite values. This need not be the case if 
the coordinate transformation that led to conformal gauge had a singularity. To cover such a case, we 
introduce a general orbifold coordinate y by 

y = 1 HKg(y)) - B , (5.11) 

where g(y) is a monotonic, continuously differentiable function for y G [2/1,2/2]- Rewriting the solu- 



tion ( |5.8P , ( |5.9[ ) in terms of this new coordinate and applying the boundary conditions (5^) leads 
to 

dsl = -K*g(y)V + K*g(yfe^ K ^ »W + g'{yfdf , (5.12) 
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where g should satisfy at the boundaries 

, n 12 , v 12 , 

S(l/l) = T^ , <7(jfe) = ±yT - (5.13) 

The upper (lower) sign applies to an increasing (decreasing) function g. Here K and tq are arbitrary 
constants [J 



Let us discuss some properties of this solution. As is well known |>^, ^3|, the metric ( 5.12J ) is 



fiat everywhere in the bulk. What makes the metric nevertheless non-trivial is the presence of the 
boundaries. While those boundaries are hypersurfaces with y = const in our coordinate frame, they 
would be mapped to de Sitter hypersurfaces in coordinates where the metric ( [5.12D takes Minkowski 
form. Indeed, if we define comoving times ti on each boundary by setting dt 2 = K 2 g(y<i) dr 2 , the 
four-dimensional boundary metrics take the form 

ds 2 ^ = -dt 2 + K 2 g( yi ) 2 e 2H * {k - Uo) dx 2 (5.14) 

with the Hubble parameters Hi given by 

H l = ±^ = ±-% . (5.15) 

12 12Mf v ; 

The above linear relation between the Hubble parameters and the potentials is quite unconventional. 
Usually the square of the Hubble parameter is proportional to the potential, as in our linear case, 
eq. (|4.6|). We should distinguish two cases for the function g. 



g(y) 7^ for all y £ [0,R]: In this case, the metric ( |5.12| ) is regular everywhere on the orbifold, as 
in our original form (|5.8| ), ( |5.9[) of the solution in conformal gauge. Since g is continuous, g{y\) 
and 5(2/2) have to have the same sign. From eq. ( |5.13 ), it then follows that 



VxVi < , Vi + V 2 < . (5.16) 

This are indeed the relations we found in conformal gauge above. 



g(yo) = for some yo £ [0, R]: In this case, the metric (15.12 ) has a horizon at y = yo- Now g(y\ 



has to be negative and 5(2/2) positive (or vice versa) and we conclude from ( |5.13| ) that 

Vi > , V 2 > . (5.17) 

This case could not be obtained from the solution in conformal gauge. 



7 Choosing g(y) = -4- — y, we obtain a form similar to the four-dimensional domain wall solution of [163, p3] and its 
five-dimensional counterpart in fed]. 
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In fig. [I] two-dimensional Minkowski space corresponding to the (r, y) plane is depicted. We have in- 
dicated the portions of this space that correspond to the above two types of solutions. As explained 
earlier, if the metric (5.12) is written in Minkowsi form by applying an appropriate coordinate transfor- 
mation the y = const orbifold planes are mapped into de Sitter hypersurfaces. In the two-dimensional 
pictures they appear as hyperbolas. The space between those hyperbolas in fig. |l] represents the orbifold 
and the lines indicate the locations of constant time, r = const. In case 1 (left figure [[]), both bondaries 
are on the same side of the light cone. Signals that travel in the bulk will always reach the boundary 
after a finite time. In particular, a signal sent from one plane will always reach the other one in a finite 
time. These causal properties are somewhat counterintuitive in that one would expect the existence of 
signals that travel exclusively in the bulk without ever hitting a boundary. The figure shows that, in 
fact, such signals do not exist. On the other hand, in case 2 (right figure [y), a signal emitted from one of 
the boundaries will never reach the other one. In this sense the two boundaries are causally decoupled. 
Again, this property is somewhat unexpected intuitively. 




y = R 



Figure 1: Parts of Minkowski space corresponding to the solution in case 1 (left figure) and case 2 (right 
figure). 

The physical length of the orbifold (obtained by integrating yjg^ over the coordinate interval [0, R], 
as usual) is static and in both cases given by 



1 1 

Rnhvs - 12 77 + 77 



i-phys 



Vi v 2 



(5.18) 
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This is positive due to the conditions fl5.16|) and ( 5.17| ), as it should. In view of eq. ( 5.15| ), this leads to 



Hi ~ — !- (5.19) 

-ftphys 

for potentials Vi of the some order of magnitude. This relation is quite promising, since it directly 
relates the Hubble parameters to the size of the additional dimension. We would like to point out that, 
unlike in the linear case, we did not assume the non-perturbative bulk potential V to depend on the 
orbifold modulus. Nevertheless, the orbifold size -R p h ys turns out to be time-independent. Moreover, it 



is fixed by eq. fl5,18|) in terms of the boundary potentials. 

We have already mentioned that the above solution is the only separating solution compatible with 
our initial assumptions of a stabilized modulus 4> and slow roll of the boundary fields fa. There is 
yet another sense in which this solution is unique. Suppose one is interested in solutions where the 
bulk moduli fields fi and eft are constant in time or slowly moving (with a negligible contribution from 
the bulk potential V to the vacuum energy), the boundary fields are slowly rolling in their potentials 
Vi and the Hubble parameter H = a changes only slowly in time. Practically, one can then neglect 
terms containing 0, (f>, fa and H in the equations of motion. Heuristically, these properties are what 
one expects from an inflating solution in five dimensions, in analogy with the four-dimensional case. 
Then, one can show that all solutions with these properties are approximated (in the sense that slow-roll 
corrections have been neglected) by ( |5 . 8| ) ( |57TO ) . To do this, one does not need the technical assumption 



of separability that we have used so far. In this sense, we have found the unique solution with boundary 
inflation in five dimensions. 

Can this solution, then, be used as the basis for an inflationary model in five dimensions? We have 
to keep in mind that we have not solved the equations of motion in general yet, but rather found a 
specific solution by imposing separability or, equivalently, "reasonable" physical conditions of what an 
inflationary solution should look like. Therefore, our solution might be very special in the sense that it, 
perhaps, can only be obtained from a set of initial conditions with measure zero. In other words, the 
solution could be unstable against small perturbations. We will analyze this question in the following 
section. 

6 Non— linear case: General solution 

In this section, we present the general solution. We start with the same setup as in the beginning of 
section 5. Since we are in the non-linear regime 

hi » i (6.i) 

we can neglect the domain wall potentials. Furthermore, we assume that <p has been stabilized by the 
bulk potential F at a point with vanishing potential energy V(4>) = (also, as a technical assumption, 

25 



we have to require the boundary potentials to be independent on (f>). Furthermore, the boundary 
potentials V% should lead to slow roll of the fields (pi. We will, however, not assume separability or slow 
time evolution of any other field. 

It is useful to introduce light-cone coordinates 

x ± = r±y (6.2) 

and rewrite the equations of motion ( [2.20| )-( |2.23] ) in terms of these coordinates choosing conformal 
gauge v = 3. Using the simplifications that follow from the above setup, one finds 

d\a - 2d+ad+(3 + d + a 2 = (6.3) 

d 2 La-2d-ad-(3 + d-a 2 = (6.4) 

d+d-a + 3d + ad-a = (6.5) 

2d + d-a + d + d-/3 + 3<9+a<9_a = . (6.6) 

The boundary conditions ( p. 26 ) and ( 2.27] ) specialize to 

e-v U = '-V U = ^ ' (6 - 7) 

where the upper (lower) sign applies to the boundary at y = y± (y = ^2)- The equations of motion ( |(0| ) 
( |6.6p are quite similar (although not identical) to those of two-dimensional dilaton gravity [70| ] with 
vanishing cosmological constant. In fact, fl6.3j)- (|6.6j ) can be obtained from the two-dimensional action 

S 2 = ~ f V^e 3a [R 2 - 24 d a ad a a] (6.8) 

with the two-dimensional metric in conformal gauge given by g a b = e 2 ^r/ a fc. Here we have used indices 
a, b, ■ ■ ■ = 0,5 for the space (r, y). Using the methods of ref. [J70|], we can find the following general 
solution of (p^|)-(p|) 

a = llnu , 3 = w — llnu , 

3 _ _ ( 6 - 9 ) 

u = u + (x + ) + n_(x ) , w = w+ (x + ) + u>_ (x ) 

where u and w are free fields, as indicated. The "left- and right-movers" u±(x ± ) and w±(x ± ) are 
not completely independent but, rather, subject to certain relations that originate from the constraint 
equations ( |6.3j ) and ( |6.4[) . Three cases can be distinguished 

• Case 1: If d+it+ 7^ and <9-_-u_ 7^ then 

w± = \Hd±u±) + C±. (6.10) 

• Case 2: If 9+u+ 7^ and cL-u_ = then 

w+ = l -Hd +u+ ) + c + . (6.11) 

and W- is arbitrary. 
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Case 3: If <9 4 



and d-U- ^ then 



W- = — In(d-u) + C- 



(6.12) 



and w + is arbitrary. 



Here C± are arbitrary integration constants. In each of the three cases, we still have two arbitrary 
functions at our disposal. While this constitutes the most general solution of the equations of mo- 
tion fl6.3| )-( |6l3|) , we have not yet taken into account the boundary conditions (|6.7|). As we will see, this 
determines one of the functions and imposes a periodicity constraint on the other one. Before we come 
to that, we observe that in case 1 we have the relation 



j3 = 2a + - ln(9<9+a<9_a) + C+ + C- 



(6.13) 



between the scale factors a and (3. Suppose that we have found a solution with an inflating three- 
dimensional universe in this case. Neglecting y-dependence, the scale factor a is then roughly given 
by a ~ Ht where H is the Hubble parameter. In this case, the second term on the right hand side of 



eq. ( 6.13 ) is approximately constant since d±a should be related to the Hubble parameter. Hence, up to 
an overall normalization we have j3 ~ 2a. This shows that, without assuming any initial fine-tuning, at 
the end of inflation, the orbifold orbifold has expanded twice as much as the three-dimensional universe. 
For the reasons discussed at the end of section 3 we, therefore, disregard this possibility. We remark 
that solving the boundary condition ( |6.7D for case 1 leads to periodicity constraints in terms of elliptic 
functions. 

Let us now turn to cases 2 and 3. Fortunately, the boundary conditions can be explicitly solved in 
these cases. We find that a and (5 can be expressed in terms of a single real function / as 

1 



a 







3 

In 



ln/(x±) 

4(/'(x ± )) 1 /2(/'( a; =F))l/2 



(6.14) 
(6.15) 



_(^l)(/(^)) 2 / 3 (/(x±))V3_ 

where /' denotes the derivative of / and, as usual, x^ = r ± y. The function / should have the 
periodicity property 



fix + 2R) 



jplRK 



(/(*))• 



-1/3 



+ k 



for all x, where 



* = !]*(-* 



(6.16) 



(6.17) 



and k is a constant. The condition ( 6.16J ) can be solved in terms of a periodic function p(x) satisfying 

p(x + 2R) = p(x) (6.18) 
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for all x. One finds 



f(x) = [p(x)e TKx + k] for if/ 

f(x) = [kx + p(x)r 3 forK = 



(6.19) 
(6.20) 



where k is another constant. The upper sign in the above solutions always corresponds to case 2, the 
lower one to case 3. The definition of K, eq. ( |6.17| ), shows that the boundary potentials need to have 
opposite sign for the solution to exist. We are, therefore, in the first case ( |5.16| ) of the previous section 
which is as expected since we have used conformal gauge. We will stick to this case, for simplicity, in 
the following. The analog of the second case ( |5.17 ) can again be obtained by employing a more general 
coordinate system. Our main conclusions apply to this case as well. Furthermore, the periodic function 
p(x) and the constant k have to be chosen such that the logarithms in ( |6.14| ) and ( 6.15| ) are well-defined. 
Apart from those restrictions, p(x) and k are arbitrary. Eq. ( |6.14 )-( |6.20D is the most general solution 



of the system (6.3)-(6.7) for the cases 2 and 3 which, as we have seen, are the interesting ones in the 
present context. Since we can more or less freely choose one periodic function we have, in fact, found a 
very large class of solutions. 

Depending on the value of K we have two different possible forms of the function /, given in 
eq. ( |6.19D and (6.20). The second option, K = 0, is realized for a vanishing "total potential energy", 
V\ + V2 = 0. Not surprisingly, in this case, the scale factor a does not inflate but (roughly) shows 
a power law behaviour, as can be seen by inserting (|6.20|) into eq. (6.14). Consequently, this second 



option is only of limited interest to us and we will concentrate on the first one, if /0. This case should 
contain our simple separating solution ( |5.8[ )-( ^10|) of the previous section. Indeed, if we choose 



p(x) = const 



k = 



(6.21) 



in eq. (|6.19| ) we recover this solution. What happens for other choices? Let us consider the case k j^ 
and p(x) periodic but otherwise arbitrary. If the argument of the exponential in eq. ( |6.19| ) is negative 
(the case that would lead to an inflationary expansion if k were zero) then, after a very short time / will 
be approximately constant, unless k is exponentially small. Inserting such an approximately constant 
/ into the expressions ( |6,14| ), ( |6.15 ) for the scale factors a and shows that the three-dimensional 
universe becomes static while the orbifold collapses. In the opposite case, where the argument of the 
exponential is positive, / falls exponentially and, hence, from eq. ( |6.14| ), the three-dimensional universe 
collapses. Either way, to have a stable configuration, we should choose k = 0. In this case, the scale 
factors read explicitly 



o 



P 



±Kx d 



Ky + ln 



hip{x~ 



t) (- 



P 1 ±- 
p 



1/2 



(±K 



-(z*) 



1/2 {p(x^)) 1 ' 2 

(p( x ±))l/2 



(6.22) 
(6.23) 
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We still have the freedom to choose the periodic function p. In the above expressions this function always 
depends on x = t ±y. It, therefore, leads to an oscillation in time with a period 2R. If we do not 
want the scale factor a to be significantly effected by this, we should choose the maximum amplitude of 
p to be sufficiently small, so that p(x) ~ const. This basically brings us back to our separating solution 
of the previous section. In essence, this solution is the only one with the desired properties within our 
setup. 

The discussion of this section has revealed two serious problems with the separating solution. First 
of all, it corresponds to a very specific choice of initial conditions satisfying either k = or k being 
exponentially small. All other values of k lead to a collapsing solution. This implies that the case k = 0, 
as it stands, corresponds to an unstable situation. A small perturbation that leads to a non-vanishing 
k will cause a collapse of the universe. The other problem is related to the presence of the periodic 
function p. This function, in fact, encodes the information about the initial inhomogeneity in the orbifold 
direction and this inhomogeneity survives the whole period of inflation. Of course, this is related to 
the fact that we are not inflating the orbifold as well which would dilute those inhomogeneities. In 
the effective four-dimensional linear case of section 4, oscillations of Kaluza-Klein modes were damped 
away quickly due to the inflationary expansion. Apparently this is no longer true in the non-linear 
five-dimensional regime. Although those modes do not affect the homogeneity of the three-dimensional 
universe directly, they could potentially be harmful. For example, eq. ( 6.22| ) shows that the "Hubble 
parameter" a contains the periodic function p{r ± y) and consequently oscillates in times. Hence, the 
modes could have some influence on density fluctuations. Also, their eventual decay into gravitons could 
leave unwanted relics. In any case, the presence of the function p contradicts somewhat the philosophy 
of inflation which is supposed to wipe out any initial information. 

Is there a possible cure for the stability problem? So far, we did not attempt to stabilize the orbifold 
in any way. It is clear that this fact is related to the instability that we encounter. While stabilizing 
the orbifold modulus in a four-dimensional effective theory by simply inventing a potential is relatively 
straightforward (although understanding the origin of this potential is not), this is not the case in five 
dimensions. The orbifold size in five dimensions is basically measured by the component 355 of the 
metric. A bulk potential for this component would break general five-dimensional covariance in a very 
strong way. Also, since 555 typically varies across the orbifold it could not be everywhere in a minimum 
of that potential. Alternatively, one could postulate a bulk potential that only depends on the zero 
mode of 555; that is, the length of the orbifold. This would break general covariance not more seriously 
than it already is broken by the presence of the orbifold. Suppose such a potential had a minimum for 
a certain orbifold length Rq. Could an orbifold modulus stabilized at this minimum be made consistent 
with our solution? The apparent problem is that the orbifold size is already fixed by eq. ( 5.1§| ) in terms 
of the boundary potentials VI. There is no obvious reason why Rq and -Rphys should, a priori, coincide. 
If we assume they do, for some reason, at a certain time, this situation could only be maintained if the 
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boundary inflatons slow-roll in a very specific way so as to leave 

R ° =u {vk + w^) (6 - 24) 

unchanged. This would correspond to a strong correlation of the motion of the two inflatons. Without 
a detailed analysis of the dynamics, which probably has to be carried out numerically, it is hard to tell 
whether this would actually happen or whether, instead, the orbifold modulus would start to strongly 
oscillate around its minimum thereby destroying inflation. 

Another way to stabilize the orbifold which avoids these problems is to have a potential for 355 on 
the boundary. Although such a potential cannot appear directly in the boundary actions in ( |2.1D , it may 
appear in one of the Bianchi identities of heterotic M-theory in five dimensions p7|, In], p3| which contain 
sources located on the orbifold plane. Particularly, a potential from gaugino condensation would mani- 
fest itself in the Bianchi identity. Generically, however, it seems to be difficult to stabilize moduli with 



potentials from gaugino condensation, in particular in the context of cosmology |34j (see however [72|). 
In ref. |7^] stabilization of the orbifold was achieved by a combination of gaugino condensation and other 
non-perturbative effects resulting from internally wrapped membranes. Although worth investigating 
in our context, all those options go beyond our simple toy model and will not be explicitly considered 
here. 

To summarize, while there seem to be interesting solutions with boundary inflation in the non- 
linear, five-dimensional regime, a closer investigation shows that they have problems related to the 
stabilization of the orbifold and to inhomogeneities in the orbifold direction that are not diluted. The 
stabilization problem is, of course, very general and we should not be surprised to encounter it in our 
cosmological context. It is conceivable that whatever eventually stabilizes the orbifold, also saves our 
non-linear cosmological solution. While the persistence of the orbifold inhomogeneities appears to be 
in contradiction with the inflationary paradigm they might, under certain conditions, be acceptable 
and give rise to interesting predictions. It certainly needs more work to finally decide these questions. 
Probably one also has to go beyond the simple model we have used in this paper. 
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